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ABSTRACT 

Context. High-mass star formation has been a very active field over the last decade, however, most studies targeted regions of lumi- 
nosities between 10"^ and 10^ L©. In contrast to that, the highest mass stars reside in clusters exceeding 10^ or even 10^ L©. 
Aims. We want to study the physical conditions associated with the formation of the highest mass stars. 

Methods. Therefore, we selected the W31 star-forming complex with a total luminosity of ~ 6 x 10^ L© (comprised of at least two 
sub-regions) for a multi-wavelength spectral line and continuum study covering wavelengths from the near- and mid-infrared via 
(sub)mm wavelength observations to radio data in the cm regime. 

Results. While the overall structure of the multi- wavelength continuum data resembles each other well, there are several intriguing 
diff'erences. The 24 /um emission stemming largely from small dust grains follows tightly the spatial structure of the cm emission trac- 
ing the ionized free-free emission. Hence warm dust resides in regions that are spatially associated with the ionized hot gas (~ 10"^ K) 
of the Hii regions. Furthermore, we find several evolutionary stages within the same complexes, ranging from infrared-observable clus- 
ters, via deeply embedded regions associated with active star formation traced by 24 yum and cm emission, to at least one high-mass 
gas clump devoid of any such signature. The ^^C0(2-l) and C^^0(2-l) spectral line observations reveal a large kinematic breadth in 
the entire region with a total velocity range of approximately 90 km s~^ . Kinematic and turbulent structures are set into context. While 
the average virial mass ratio for W31 is close to unity, the line width analysis indicates large-scale evolutionary diff'erences between 
the southern and northern sub-regions (GIO. 2-0.3 and GlO.3-0.1) of the whole W31 complex. A color-color analysis of the IRAC 
data also shows that the class II sources are broadly distributed throughout the entire complex whereas the Class O/I sources are more 
tightly associated with the active high-mass star-forming regions. The clump mass function - tracing cluster scales and not scales of 
individual stars - derived from the 875 jum continuum data has a slope of 1.5 + 0.3, consistent with previous cloud mass functions. 
Conclusions. The highest mass and luminous stars form in highly structured and complex regions with multiple events of star forma- 
tion not always occurring simultaneously but in a sequential fashion. Warm dust and ionized gas can spatially coexist, and high-mass 
starless cores with low-turbulence gas components can reside in the direct neighborhood of active star-forming clumps with broad 
line- widths. 
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1. Introduction 

High-Mass star formation research has made significant progress 
over the last decade, from a theoretical as well as an observa- 
tional point of view. For recent reviews covering various as- 
pects of high-mass sta r formation, we refer to, e.g ., Beuthe r et al. 
(|2Q07l);|Ziiinecker & Yorke (2007l);[Bonnell et al.K2007 ) ; McKee 



|& Qstriker|p007| ); |Krumholz & Bonnell|p007| ). Unfortunately, 



because of statistical selection eflfects and general low-number 
statistics at the high-luminosity end of the cluster distribution, 
one observational drawback so far was that most (sub)mm stud- 



< 10^ Lo, hence stars of mostly less than 30 M© (e.g.,|Molinari 


et al."1996'i "Plume et al.pi997 


, ISridharan et al.||2002[ |Mueller 


et al. 2002; Fontani et al.pOOS 


). However, some notable excep- 



tions exist, e g., GlO.6-0.4 (Keto, 2002), G3 141+0.31 ( [Cesaroni 
leTaLl [1954^ or W49 ffiomeier & Alvesj [2005] ). Hence, we 
are lacking solid observational constraints on the physical and 
chemical conditions of young high-mass star-forming regions 
forming the most luminous and high-mass stars within our 
Galaxy. In an eflfort to overcome these limitations, we selected 
one of the most luminous (~ 6x10^ L©) but not too distant 



(~6kpc compared to, e.g., ~11.4kpc for W49) Giant Molecular 
Cloud/high-mass star formation complex G10.2/G10.3 (also 
known as W31) for a detailed multi- wavelength investigation 
based on APEX ^^CO and C^^O mapp ing observations, th e 
875 yum ATLASGAL continuum survey ( Schuller et al.|[2QQ9l), 
Spitzer GLIMPSE/MIPS G AL mid-infrared data ( [ChurchwelT 



et al.| 2009 [ Carey et al. 2009), and previously published cm 
continuum data ( |Kim & Koo||2002| ). 

The high-mass star formation complex W31 contains the 
sub-regions GlO.2-0.3 and GIO. 3-0.1 which have already been 
identified as Hii regions and IRAS sources (IRAS 18064-2020 
and IRAS 18060-2005) with larg e luminosities (e.g., . Wood ward] 
let al.||1984t [Ghosh et all|1989| ). Its dist ance is much debated 
and literature values vary between 6 k pc (|Wilson[ 1974[|Downes 
eFaL 1980), 14.5 kpc ( [Corbel et al.||T997 |) and 3.4 kpc (|Blum 
et al.[ [200 1| ). The two main regions are separated by approxi- 
mately 15' (Fig. 1). At an assumed distance of 6 kpc, the total 
luminosity and gas mass of the regio n amounts to ~ 6 x 10^ L© 
and ~ 6 X 10^ M© ( Kim & Kool |2002 ), respectively. The individ- 
ual luminosities of the main far-infrared peaks of the two sub- 
regions GlO.2-0.3 andGlO.3-0.1 a re estimated to be ~ 8 x 10^ L© 
and ~ 6 X 10^ L©, respectively (iGhosh et al.l |l989|). Within 
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the larger-scale Hii regions, two ultracompact Hii regions were 
identified GlO.15-0.34 and GlO.30-0.15 (Wood & Churchwell, 
^1989 ). The luminosities of these two UCHiis based on IRAS 
dataTare ~ 1.5 x 10^ Lp and ~ 7 x 10^ L©, respectively ( Wood 
|& Churchwell[ T989]), con sistent within a factor 2 with the es- 
timates from [Ghosh et al.| ( |1989| ). It is interesting to note that 
the infrared clusters in both regions (GlO.2-0.3 and GIO. 3-0.1) 
are spatially off'set from the UCHii regions, indicating diff'er- 
ent episodes of high-mass star formation. Furthermore, Walsh] 
et al. (1998]) identified four distinct Class ii CH3OH maser po- 



continuum obser vations with a spatia l resolution of ~ 31" x 25' 
first pubHshed by |Kim & Kool ( [j502l ). 



sitions toward the northern GlO.3-0.1 sub-region also known as 
IRAS 18060-2005) but none toward the southern GlO.2-0.3 sub- 
region (a.k.a. IRAS 18064-2020). 

While small submm continuum maps of individual clumps 
within each of the regions had been obtained over the last few 
years (e.g., the SCAMPS project, [ Thomp son et al.|2006|), there 
did not exist (sub)mm continuum data encompassing the whole 
GMC complex. With the advent of the 875 //m survey of the 
southern Galactic plane (ATLASGAL, |Schuller et al.|2009| ), we 
now have a complete census of the submm continuum emis- 
sion of this large-scale star-forming region. Kim & Koo, (2002) 
mapped the whole region in ^^CO(l-O) and CS(2-1) at a rela- 
tively low spatial resolution of 60''. They reveal many interesting 
large-scale features (e.g., that G10.2 and GIO. 3 may be located 
on a shell-like structure). Similarly, [Zhang et aL] ( |2005| ) observed 
the complex in C^^O(l-O) with a poor grid separation of T. All 
these studies are not resolving the small-scale structure, and they 
are not sensitive to the higher-density gas components. 



2. Observations and archival data 

The Ci^0(2-1) and ^^C0{2-\) data at 219.560 GHz and 
220.399 GHz were observed simultaneously with the Atacama 
Pathfinder Experiment (APE^Q Giisten et al. 2006) in July 
2009 in the 1 mm band over the entire W3 1 complex in the on- 
the-fly mode. The APEXl receiver of the SHeFI receiver fam- 
ily has receiver temperatures of -130 K at the given frequency 
( [Vassilev et a l. , 2008 ), and the average system temperatures dur- 
ing the observations were ~250 K. Two Fast-Fourier-Transform- 
Spectrometer (EFTS, Klein et al. 2006| ) were connected covering 
~2GHz bandwidth between 219 and 221 GHz with a spectral 
resolution of ~0.17 km s"^ The data were resampled to 1km s"^ 
spectral resolution and converted to main beam brightness tem- 
peratures with forward and beam efficiencies at 220 GHz of 
0.97 and 0.82, respectively ( [Vassilev et al. 2008[ ). The average 
Icr rms of the final spectra is -0.98 K in Tmb. The OFF-positions 
is approximately 42724' offset from the center of the field (se- 
lected via t he Stony Brook Galactic Plane CO survey, Sanders 
et al.|1986] ), and it apparently has emission in a velocity regime 
between 27.5 and 36.5 km s"^ showing up as artificial absorption 
features in our data (see section |T2| ). The FWHM of APEX at 
the given frequencies is ~ 27.5". 

The 875 yum data for the region are from the APEX 



ATLASGAL survey of the Galactic plane ( [Schuller et aL|[20Q9] ). 
The Icr rms of the data is ~70mJybeam and the FWHM ~ 
19.2". The Spitzer IRAC 3.6, 4.5, 5.8 and 8.0yum images as weU 
as the MIPS 24 jim are taken from the GLIMPSE and MIPSGAL 
Galactic plane su rveys respectively ( [Churchwell et al.j |2009[ 
Carey et aL) 2009]). Furthermore, we employ the 21 cm radio 



3. Results 

3.1. Multi-wavelength continuum Imaging 
3.1 .1 . General properties 

Figure 1 shows a compilation of the different continuum datasets 
analyzed in this work ranging from the mid-infrared with IRAC 
8 yum and 24 [im data, to the submm regime at 875 [im and then 
to radio wavelengths at 21 cm. In a simplified picture, we se- 
lected the mid-infrared images to trace warm dust emission, the 
submm data to study cold dust emission, and the cm observations 
to investigate the ionized gas components. We clearly identify 
the large-scale structure of the two star formation and Hii region 
complexes GIO. 2-0.3 in the south and GlO.3-0.1 in the north. 
At all wavelengths the two complexes exhibit significant sub- 
structure, probably most strongly recognizable in the 875 yum 
image. Since we are dealing with two very luminous but already 
relatively evolved star formation complexes it is not surprising 
that many emission features are visible at all presented wave- 
lengths. They show ionized gas as well a s warm dust emission 
caused by the luminous O and B stars (Blu m et aL) |2001[ |Bik 



|et aH |2005|), and furthermore a large amount of cold dust emis- 
sion from the still existing gas/dust reservoir present in the whole 
region. While the former signifies that star formation is already 
ongoing for a considerable amount of time (> 10^ - 10^ yr), the 
latter shows that the original gas cloud is not yet dispersed and 
may still allow further accretion. 

Furthermore, the 24yum and 21 cm data spatially resemble 
each other well, for example the cometary shape of the north- 
ern region GlO.3-0.1. An even clearer representation of the cor- 
relation between 21 cm and 24yum emission can be achieved if 
one plots the observed fluxes pixel by pixel. Figure |2] presents 
such a "scatter plot" correlating all pixels with 24 yum flux above 
120MJy sr"^ (corresponding to the approximate outskirts of the 
Hii regions) and the 21 cm emission above the 3cr level of 
-O.OlJybeam"^ Only scales of the order the beam size of the 
21 cm observations were correlated (2-pixel steps correspond- 
ing to 29"). This figure outlines the good spatial correlation be- 
tween warm dust and ionized gas emission over a broad range 
of fluxes in the W3 1 complex. Hence warm dust below the sub- 
limation temperature of ~1500K spatially coexists wi th much 



^ This publication is based on data acquired with the Atacama 
Pathfinder Experiment (APEX). APEX is a collaboration between 
the Max-Planck-Institut fur Radioastronomie, the European Southern 
Observatory, and the Onsala Space Observatory. 



warmer ion ized gas of temperatures around 10 " ^ K (e. g., |Helfand 
eraLl[20Q6| ). Recently, [Everett & Churchwell] ( [2QT0| ) argue that 
such dust may stem from continuous ablation processes of small 
cloudlets that have not been entirely destroyed by the Hii region 
yet. However, the spatial resolution of the 21 cm observations 
corresponds to linear scales of about ~ 0.9 pc. At these scales it 
is also possible that the observed 21 cm/24 jim correlation stems 
from an interface between the edges of Hii regions and dense 
dust shells which are just not resolved by the observations. One 
should keep in mind that the 24 yum band is mainly sensitive to 
very small, s tochastically heated dust grains (e.g., Draine|2003j 
Carey et ar| |2009) and that, therefore, strong 24/jm emission 
does not necessarily imply large dust column densities. 

While in many cases we see clear associations of emission 
at difl'erent wavelengths, we also witness exactly the contrary. 
Figure |3] presents a zoom into the northern complex GlO.3-0.1, 
where we also added the near-infrared K-band data from Bik] 
( |2004| ). The latter outline the location of the infrared cluster 
that is ofl'set from the emission we see at other wavelengths. 



H. Beuther et al.: High-mass star fomation at high luminosities: W31 at > 10^ L( 



3 



O 




18'09"50' 18'09"00' 18'09"50' 18'09"00' 1 8'09"50' 18'09"00' 

R.A. (J2000.0) 



Fig. 1. Continuum images of the W31 region. The shown wavelengths are marked in each panel. The 875 yum data are contoured 
in 3cr steps from 3 to 12cr and continue in IScr steps from 15cr onwards (Icr ^ 70mJy beam"^). The 21 cm image is contoured 
in 3(T steps from 3 to I2cr (Icr ^ SmJybeam"^), continuing in 0.24 Jybeam"^ steps from 0.12 to 0.84beam"\ and then go on in 
0.48 beam"^ ste ps. The 24 //m m ap is saturated toward the peak positions. In the top-left panel the 3-pointed stars mark the po sition 
of an 0-star by ( jBik et al. |2005 ) and the approximate center o f the southern infrared-cluste r discussed by Blum et al. (2001 ). The 



two white 5 -pointed stars show the positions of UCHii regions ([Wood & Churchwell 



velocity-wise is associated with cloud complexes outside our field of view (Sec. 3.2). A scale-bar is shown in the top-left panel. 



1989|). The ellipse marks the emission which 



The 875 yum and 21 cm beam sizes are shown in the top-right corners of the top-left and bottom-left panels, respectively. In the 
bottom-left panel, the triangles and stars show IRAC-identified protostellar class O/I (the two classes are combined) and class II 
candidates, respectively, following [Allen et al.| ( [2QQ4] ); |Megeath et al.| ( [2Q04| ); |Qiu et aTlp008| ). 
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Fig. 2. Correlation between 24yum fluxes and 21cm fluxes 
for the W31 complex. Only pixels with lAjim flux above 
120MJysr"^ (corresponding to the approximate outskirts of 
the Hii regions) and 21 cm emission above the Scr level of 
~0.01Jy beam"^ are plotted. Furthermore, only correlations on 
scales of the beam size of the 21 cm observation were used (2- 
pixel steps corresponding to 29'')- 



The 0-star marked in Fig. [T] (| Bik et alT) |2QQ5| ) is part of that 
cluster. Directly east of the infrared cluster we find the associ- 
ated submm clump G10.3E (clump 4 in Table [2]) and further to 
the north G10.3NE (clump 2 in Table |2]), both associated with 
mid-infrared, cm and Class ii CH3OH maser emission indicat- 
ing active star formation. Moving west from the infrared cluster 
we identify a clear emission peak at all wavelengths (G10.3C 
in Fig. [3] and clump 1 in Table [2]) marking already a relatively 
evolved part of the UCHii region. Going west, the next source 
G10.3W1 is still associated with mid-infrared, cm continuum 
and Class 11 CH3OH maser emission (G10.3W1 and G10.3W2 
merge into clump 6 in the smoothed data of Table [2]). Although 
at a weaker level, this also indicates active high-mass star for- 
mation from this gas clump. However, in strong contrast to these 
sub-regions, the most western submm peak G10.3W2 in Figure 
|3] shows comparably strong submm continuum emission but a 
complete absence of any 24 jim and 21 cm emission. Hence this 
is a high-mass gas clump at a very early evolutionary stage, po- 
tentially still in a starless phase prior to any star formation activ- 
ity. The projected separation of the G10.3W1 and G10.3W2 peak 
position of ~ 42'' corresponds at the assumed distance of 6 kpc 
still to a projected spatial linear separation of -1.2 pc. Hence 
each of the submm emission peaks covers scales correspond- 
ing to clusters, and therefore each may form a smaller cluster - 
likely also containing high-mass stars - within the large-scale 
environment of W3 1 . 

Figure [T] also shows the distribution of protostellar class O/I 
(the two classes are combined) and pre-main- sequence class II 
sources identified by Spitz er IRA C data based o n the selection 
crit eria develop ed by Alle n et al.| pQQ4), Megea th°eFaL| ( |20Q4| ) 
and Qiu et al.] ( 2Q08j . While the non-detection of any source 
toward the centers of the two sub-regions is likely an artifact 
due to saturation and confusion in these areas, one can tenta- 
tively identify a trend that the class O/I sources are more closely 
associated with the centers of activity than the class II sources 
which appear more widely distributed over the field of view. 



Quantitatively speaking, the mean separations of the northern 
class O/I and class II sources from the UCHii region G10.30- 
0.15 are ~ 321'' and ~ 428", respectively, whereas the corre- 
sponding mean separations for the southern sources with respect 
to the UCHii region GlO.15-0.34 are ~ 936" and ~ 1048", re- 
spectively. These data also indicate that the southern cluster is 
spatially more distributed and may potentially be already more 
evolved (see section [4TT]). 



3.1 .2. Clump properties of the 875 yum continuum emission 

We can now locate the dense 875 jim gas and dust clumps at 
a spatial resolution of 19", corresponding to linear scales of 
-0.5 pc at an assumed distance of 6 kpc. While two of the dense 
cores are associated with the ultracompact Hii regions G10.15- 
0.34 and GlO.30-0.15 (Fig.[T]), there are obviously many dense 
gas clumps that are promising candidates of ongoing and poten- 
tially even younger high-mass star formation. At the given spa- 
tial scales of our resolution limit, each submm sub-source does 
not form individual stars but they all should be capable to form 
sub-clusters within the larger- scale region. 

To systematically identify gas and dust clumps and to study 
their properties in such a large and complex reg ion we employ 
the clumpfind algorithm ( Willia ms et al.| |1994| ) with 3(T con- 
tour levels (210 mJy beam"^). With this procedure we identify 73 
submm continuum clumps over the entire W31 region (Fig. [4]). 
These clumps are listed in Table[T]with the spatial coordinates of 
their peak positions, their integrated and peak fluxes, and their 
linear efl'ective radii reff. Assuming that the submm emission 
stems from optically thin dust, we can calculate the H2 gas col- 
umn densities and masses following ,Hildebrand| (| 1 9 8 3 1) ; |B euther 



et al.| ( [20Q2||2005| ). Since we do not know the temperature sub- 



structure of the entire complex, we use an average dust tempera- 
ture of 50 K that should be a reasonable proxy for such an active 
region (e.g., Sridharan et al.|2002| . Furthermore, we use a dust 
grain emissivity index ofp = 2 (corresponding to a dust absorp- 
tion coefficient of /cgvs ^ 0.8 cm^ g~^), and a gas-to-dust m ass ra- 
tio of 186 is assumed f Jenkins| [2504l [Draine et"aLl [2007 ). Table 
[T] lists also the derived H2 column densities and masses of the 
respective gas clumps. With the given uncertainties of tempera- 
ture and dust composition, we estimate the masses and column 
densities to be accurate within a factor ~3. 

The derived clump masses range between 120 and 8200 M© 
with a total gas mass in the region of approximately 1.2x10^ M©. 
The peak column densities vary between 1.5 x 10^^ and 3.4 x 
10^^ cm"^. While th e proposed threshold for hig h-mass star for- 



2008b corresponds to 



mation of 1 g cm"^ ( [Krumholz & McKee 
column densities of ~ 3 x 10^"^ cm"^, this does not imply that 
most clumps are not capable of high-mass star formation, the 
data clearly show the opposite. In contrast to this, the calcu- 
lated column densities are measured with a linear beam size of 
-0.5 pc, and hence only the average column densities over such 
area have the derived values. The intrinsic column densities at 
smaller spatial scales are significantl y higher assuming a ty pi- 
cal density distribution oc r"^ (see also Vasyunina et al.|20Q9 for 
similar estimates conducted for infrared dark clouds). 

This implies that a large number of the gas clumps in this re- 
gion should be capable of forming sub-clusters containing high- 
mass stars. To estimate how much mass a gas clump needs to 
have to form a high-mass star, we produce stellar cluster toy- 
models assuming a star formation efficiency of 30% and an ini- 
tial mass function (IMF) following Kroupa, ( ,2001) . To form at 
least one high-mass star of, e.g., 20 or 40 M©, in this scenario the 
initial gas clumps have to have masses of ~ 850 or ~ 1900 M©, 
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Fig. 3. Zoom into the G10.3 complex. In all 4 panels the solid contours present the 875 yum continuum contours (start at Scr and 
continue in 9cr steps, Icr ^ TOmJybeam"^) whereas the grey-scale shows other wavelength data as outlined over each panel. The 
dashed 21 cm contours in the right panel go from 15 to 60mJybeam"^ in 15 mJy beam"^ steps, and from 120 to 840 mJy beam" ^ 
in 240 mJy beam" ^ steps. The K-band data are taken from Bik (2004). The white central part of the 24yum image is saturated. The 
triangles in the left panel mark th e Class ii CH3QH maser posi tions from |Walsh et al.||l998 ) . The white star in the right panel marks 
the position of an UCHii region ( [Wood & ChurchwelT 1989), and the additional labels there name sources discussed in the main 
text. 



respect ively. Similar results w ere recently obtained observation- 
ally by Johnston et al. (2009). This clearly indicates that a sig- 
nificant fraction of gas clumps in the W3 1 region is capable to 
form stars in excess of 20 M©. 



Using the clump masses from Table [T] we also tried to de- 
rive a clump mass function AN/ AM oc for the region. To 
overcome fitting artifacts due to diff'erent bin sizes, we fitted the 
data, systematically changing the bin width (see also Rodon et 
al. subm. to A&A). To avoid any sensitivity cutoff' we fitted only 
clumps above 1000 M© and furthermore required at least three 
non-empty bins. This procedure resulted in 36 fits to the data 
with diff'erent bin sizes each time, allowing us to asses better 
the error margins of this approach. Figure [5] shows one example 
fit. While the Poisson errors for each individual fit are relatively 
large, the combined assessment of all diff'erent fits from varying 
bin sizes in this fitting procedure as well as a comparison to a cu- 
mulative fit, allows us to derive a power-law distribution index 
within reasonable error margins. Our derived power-law index a 
from this procedure is ~ 1 .5 ± 0.3 (see discussion in section 4.2 ). 



3.1 .3. Smoothed 875 yum data 

For a better comparison with the spectral line data discussed in 
the following section, we also smoothed the 875 fim continuum 
data to the same spatial resolution as the spectral line obser- 
vations (~ 27.5'')- On this smoothed map with a lower rms of 
50mJybeam"^ we again apply the clumpfind algorithm to ex- 
tract sources, resulting in 44 clumps compared to the 73 clumps 
in the original data (see Fig.[4]for comparison). For these clumps 
we also calculated masses and column densities as conducted in 
the previous section for the un- smoothed data. In addition for 
comparison purposes with the spectral line data, we also cal- 
culate the mass just within the central beams toward each peak 
position Mpeak- All the parameters extracted for the smoothed 
clumps are listed in Table |2] 

3.2. Kinematics from spectral line data 

Star formation processes significantly shape the kinematic and 
dynamic properties of the molecular gas. While early-on dense 
starless gas and dust clumps usually exhibit narrow line widths 
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Fig. 4. The color-scales show the clump structure and boundaries derived with the automized clumpfind procedure on the original 
(left) and smoothed (right) 875 yum continuum observations. The contours show the original 875 yum observations with the same 
contour levels as in Fig. [T] The wedge reflects the clump numbers as in Tables [T] & |2] A few discussed clumps are labeled in the 
right panel (see also Fig.fTl]). 



because no (proto) stellar feedback has yet altered the origi- 
nal gas properties, during ongoing star formation, the natal gas 
clumps are strongly influenced by the central star-forming pro- 
cesses which is reflected in the observable spectral signatures. 

The ^^C0(2-l) and C^^0(2-l) observations allow a kine- 
matic analysis of the region. Figure |6] presents the spectra of 
both species averaged in each case over the entire field of view 
of the observations. In particular the ^^C0(2-l) data show the 
very broad velocity range present in the region. The figure also 
marks diff'erent velocity regimes used by us to produce moment 
maps of the data (see below). The velocity regime from 27.5 
to 36.5 km s"^ appears in absorption in this averaged spectrum. 
This, however, is an artifact because we see it in all spectra and 
hence it stems from emission at this velocity in the OFF position. 



Interestingly, this velocity range is exactly the same as the ve- 
locities we find emission for toward the submm continuum peak 
marked by an ellipse in Figure [T] Hence, this continuum feature 
in our maps should be associated with other cloud complexes 
outside of our observed field of view. 

Figure [t] presents ^^C0(2-l) integrated intensity images of 
diff'erent parts of the ^^C0(2-l) spectrum from Figure [6] We 
omit the velocity structure lower than -6kms"^ because this 
is only lower-intensity more diff'use emission which is hard to 
properly image. However, it very likely belongs to the overall 
W3 1 region. The main velocity component present in both W3 1 
complexes (GlO.2-0.3 and GlO.3-0.1) is the main peak between 
2 and 21kms"^ In contrast to this, the two velocity compo- 
nents between -6 and 2kms"^ and 36.5 to 58kms"^ are clearly 
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Fig. 7. The grey-scale shows different ^^C0(2-l) integrated intensity images where the integration regimes are shown in the top 
of each panel. Contour levels go from 20 to 80% in 20% levels of the respective peak intensities. The dotted contours outline the 
875 yum emission at a 3cr level from Fig.[T] 



only associated with specific sub-regions within G 10.2-0. 3 and 
GlO.3-0.1, respectively. As outlined in the previous section, the 
fact that we see the feature between 27.5 and 36.5 km s"^ not 
at any other velocity but only as a negative feature throughout 
the map due to emission from the OFF-position supports the 
idea that this structure is likely associated with other clouds out- 
side our field of view. Similarly, the structure toward the north 
between 58 and 77kms"^ may also not be associated with the 
W31 region but could be at a diff'erent distance. While the main 
component between 2 and 21 kms"^ is widely distributed, it is 
difl&cult to safely distinguish whether the other velocity compo- 
nents in the region are simply chance alignments of projected 
sources at diff'erent distances, or whether there exists an addi- 
tional (projected) large-scale velocity gradient from blue-shifted 
emission toward the south-east to more red-shifted emission to- 
ward the north-west. As visible in Fig.|6j the averaged C^^0(2- 
1) spectrum exhibits in general similar velocity features as that 
of ^^C0(2-l), however, at a lower level. Therefore, in Fig. [s] we 
only show the integrated C^^0(2-l) emission of the main veloc- 
ity component between 6 and 19 km s"^ 

Figures [8] and [9] present the 0th, 1st and 2nd moment maps 
(integrated intensity, intensity-weighted velocity structure and 
intensity-weighted line-width structure) of the main velocity 
components of C^^0(2-l) and ^^C0(2-l), respectively. While 



the general structure of the ^^CO and C^^O emission of this main 
velocity components agrees relatively well with the 875 yum con- 
tinuum emission, in particular the integrated ^^C0(2-l) maps 
exhibits much more extended gas structure than the dust contin- 
uum emission. This is likely a combined eff'ect of the limited sen- 
sitivity of the ATLAS GAL data on the one hand (approximate 
rms of 50mJy beam-\ [SchuUer et aL]|2009t , and on the other 



hand it may also be due to the general large-scale spatial filter- 
ing eff'ect inherent to any bolometer observations. Furthermore, 
the C^^O integrated emission spatial structure agrees well with 
the high column density part of the dust continuum maps out- 
lined by the red contours in Figure [8] 

While some sub-regions exhibit additional velocity compo- 
nents (see Fig.JTjand discussion above), the 1st moment maps 
in Figs. [8] and |9[do not exhibit a clear velocity diff'erence be- 
tween the two large-scale sub-complexes within the main ve- 
locity component. Rather in contrast, both regions - GlO.2-0.3 
and GlO.3-0.1 - show velocity sub-structure associated with 
the width of the main velocity component. Probably more im- 
portant, the 2nd moment maps show a very broad distribution 
of line-widths, extending in ^^C0(2-l) to values in excess of 
lOkms"^ Although still comparably narrow (see Av for clump 
1 in Tables [T]& [2]) in the northern region (GIO. 3-0.1), the broad- 
est line- width there is measured toward clump 1 (Table[2]), which 
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Fig. 8. The grey-scale shows C^^0(2-l) moment maps: 0th moment or integrated intensity on the left, 1st moment or intensity- 
weighted velocity in the middle and 2nd moment or intensity- weighted line-width to the right. The black contours in the left panel 
show the integrated C^^0(2-l) data and go from 20 to 80% in 20% steps. The dotted contours outline the 875yum emission at a 
3(T level from Fig.[T] and the red contours show additionally the high-intensity 875 yum emission in IScr steps from 15c r onwards 



(Icr ^ TOmJybeam ^). In the right panel, the 3-pointed stars in the north and south mark the position of an 0-star by Bik et al 
( 2005|) and the center of the southern cluster discussed by |Blum et al.| ( |2001| ). The 5-pointed stars mark the two UCHii regions 
jWood & Churchwell|[T989l ). 



is spatia lly associa t ed wit h the near-infrared cluster around the 
0-stars ( Bik et al. 2005| ), hence the most evolved part of this 
sub-complex. Although the southern region Gl 0.2-0. 3 exhibits 
in general much broader line widths compared to the northern 
region, there the broadest C^^0(2-l) and ^^C0(2-l) line-widths 
are attributed again t o regions in the close vicinity of the infrared 
cluster discussed by Blum et al. ( 2001| ). Toward the southern 
and northern regions, the UCHii regions both exhibit compara- 
bly broad line width, however, always narrower than toward the 
gas clumps associated with the infrared clusters. 



3.3. Line parameters and virial analysis 

In particular the C^^0(2-l) spectral line data are suited for 
a more detailed spectral analysis. Therefore, we extracted the 
C^^ 0(2-1) spectra toward all peak positions identified by the 



clumpfind procedure for the dust continuum data. However, one 
should keep in mind that the spatial resolution of the continuum 
and line data is diff'erent, 19.2'' and 27.5'', respectively. While 
for completeness, we extract the spectra toward all positions 
from the high- spatial-resolution continuum data (Table [T]), for a 
better comparison, it is more useful to extract the line parameters 
toward the positions extracted from the smoothed 875 yum con- 
tinuum map (Section [3.1.3| ). Table[2]lists the derived line param- 
eters toward the smoothed dust continuum peak positions from 
Gaussian fits to the C^^ 0(2-1) line profiles. Line parameters are 
peak temperature Tpeak, integrated intensity J T^Q^^dv, peak ve- 
locity Vpeak and Full Width Half Maximum line width Av. 

These line parameters can be used to derive physical quan- 
tities, in particular H2 column densities and virial masses. To 
calculate the C^^O and corresponding H2 column densities, 
we followed the standard Local Thermodynamic Equilibrium 
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Fig. 9. The grey-scale shows different ^^C0(2-l) moment maps: 0th moment or integrated intensity on the left, 1st moment or 
intensity- weighted velocity in the middle and 2nd moment or intensity- weighted line-width to the right. The black contours in the 
left panel show the integrated intensity ^^C0(2-l) data from 20 to 80% in 20% steps. The dotted contours outline the 875 /im 
emission at a 3cr level from Fig.[T] and the red contours show additionally the high-intensity 875 jim emission in I5cr steps from 
15(T onwar ds (Icr ^ 70mJy beani^). In the right panel, the 3-pointed stars in the north and south mark the position of an 0-star as 
reported by Bik et al.|(|2005|) and the center of t he southern cluster discussed by |Blum et al.| ( [2001| ). The 5-pointed stars mark the 
two UCHii regions ([Wood & ChurchwelTlp^89l). 



(LTE) analysis ( [Rohlfs & Wilson[ [20061 ) assuming optically thin 
C^^0(2- l) eni ission again at an average temperature of 50 K 



(section |3.1.2|). Comparing individual ^^CO and C^^O spec- 



tra of the region indicates that this assumption is reasonable. 
The resulting H2 column densities are listed in Tables [T] and 
|2] While higher/lower temperatures for individual sub- sources 
would lower/raise the column density estimates, we use again a 
uniform temperature to better compare later with the dust contin- 
uum derived results (Fig. 10). Non-optically thin C^^O emission 
would also raise the derived column density values. 

Furthermore, we calculated virial masses for th e given 
C^^0(2-l) line widths following [MacLaren et aL](|l988t 



(1) 



with k2 = 126 for a density profile p oc r^^, R the radius of 
the clump defined as half of the FWHM of the beam and Av the 
line width given in Table |2] Under the assumption of a flatter 
density profile (p oc r"^) or a Gaussian density distribution the 
viri al masses wou ld be app roximately factor s of 1.5 or 3 higher 
(Mac Laren et aH [l988; Si mon et al.T|2001| ). We assume a fac- 
tor 2 uncertainty for the calculated virial masses. The derived 
masses are listed also in Tables [T] and O 

Figure [T0| presents a comparison of the H2 column densi- 
ties derived on the one hand via the 875 yum continuum emission 
and on the other hand from the C^^0(2-l) emission at the same 
spatial resolution of 27.5''. Although the real errors may even 
be larger (see discussions above), for clarity we draw "only" 
25% error margins in the figure. For almost all positions, the 
H2 column densities derived from the 875 yum data exceed those 
derived via the C^^O emission. Even considering the given un- 



10 



H. Beuther et al.: High-mass star fomation at high luminosities: W31 at > 10^ L( 




Fig. 5. Example clump mass function for the W3 1 complex for 
one bin size (here logarithmic bins with 10^*^-^ < M < 10^^+^^*^-^ 
with k din integer, 36 different bin sizes have been used to over- 
come fitting artifacts). Only clumps with masses > lOOOM© 
were used for the fit. 




from aibfim (cm ) 



certainties for the diff'erent column density derivations, this is 
plausible since with a critical density of ~ 10"^ cm"^, C^^ 0(2-1) 
does not trace the highest density gas and therefore may miss the 
highest column density regions. Similar results were recently ob- 
tained by |Walsh et aL] ( [2OTQ| ) for the regions NGC6334I & I(N) 
where even the very optically thin CO isotopologue C^^O(l-O) 
did not trace the dust continuum derived column density peak of 
the region. 



A diff'erent comparison is presented in Figure 1 1 where we 
show the gas masses derived from the 875 yum dust continuum 
peak fluxes Mpeak against the masses obtained from the C^^0(2- 
1) line width under the assumption of virial equilibrium Myir. 
Similarly to the column density figure, we draw "only" 25% 
error margins for clarity reasons. To first order, we find no 
clear trend but a "scatter plot". However, excluding clumps 1, 
2, 4, 6 and 10 (Fig. [4] and Table |2]), although with a consid- 
erable scatter most other clumps are located in the vicinity of 
the Myir/Mpeak ~ 1 relation. Therefore, many of these clumps 
may not be far from virial equilibrium. Inspecting clumps 1 , 2, 
4, 6 and 10 below the 1:1 relation in Figure [TT] in more detail, 
interestingly we find that all these sources are associated with 
the northern sub-region GIO. 3-0.1 (Figs 11 4]&|3]), whereas the 
southern region GlO.2-0.3 is dominated by higher Myir/Mpeak 
ratios. Hence, there appears to be a significant diff'erence in the 



Fig. 10. Comparison of H2 column densities derived from the 
875 yum dust continuum and C^^0(2-l) at the same spatial reso- 
lution of 27.5''. We draw 25% error-bars although the real errors 
may be even larger as discussed in the main text. The solid line 
marks the 1 : 1 relation. 



turbulent line width contribution between both sub-regions (see 
discussion in section 14. lb. 



4. Discussion 

4.1. Different evolutionary stages and kinematic properties 

As already mentioned in the introduction, the spatial off'set be- 
tween near-infrared clusters and UCHii regions toward the north- 
ern GlO.3-0.1 and the southern GlO.2-0.3 complex are indica- 
tive of several episodes of high-mass star formation w ithin e ach 
of the sub-regions. Furthermore, as outlined in section 3.1.1 and 
Fig.|3j within the GIO. 3-0.1 region, we do not only have the in- 
frared cluster and the UCHii region, but we also find at least one 
high-mass gas clump G10.3W2 without any cm or mid-infrared 
emission indicative of star formation. Hence this clump may still 
be in a starless phase prior to active star formation. Therefore, 
from an evolutionary point of view, the northern GlO.3-0.1 com- 
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Peak gas masses derived from 875/2m dust data (M^^J 

Fig. 11. Comparison of masses derived from the 875 yum dust 
continuum as well as via the C^^0(2-l) line width under the as- 
sumption of virial equilibrium at the same spatial resolution of 
27.5''. We draw 25% error-bars although the real errors may be 
even larger as discussed in the main text. The solid line marks the 
1 : 1 relation. The numbers in the bottom-right part correspond to 
the clumps discussed in the main text and marked in Figure |4] 



plex hosts several high-mass star-forming gas and dust clumps 
- all potentially capable to form high-mass clusters - in at least 
three different evolutionary stages: a more evolved infrared clus- 
ter where the main sequence stars have already emerged, high- 
mass protostellar objects with and without embedded UCHii re- 
gions that may still be in an ongoing accretion phase, and high- 
mass starless clump candidates (see Table [T] for their parame- 
ters). Although there is no robust proof that the infrared cluster 
actually triggered the star formation processes in the neighbor- 
ing clumps, the spatial structure with the infrared cluster almost 
at the geometrical center of GlO.3-0.1, then the high-mass pro- 
tostellar objects following in the surrounding, and the high-mass 
starless clump candidates at the edge of the region, indicates that 
triggered star formation may be a possible cause for the different 
evolutionary stages in this region. 

Toward the southern GlO.2-0.3 complex we also identify an 
infrared cluster and several high-mass protostellar objects, but 
there we do not as unambiguously find regions classifying as 
high-mass starless clump candidates. For a broader discussion 
of evolutionary s equences in high-mass star formation, we refer 
to the reviews by [Beuther et al.| pOQ7| ) and |Zinnecker & Yorke 
( |2QQ7l ). 

In addition to the evolutionary differences within each of the 

lis 



two main regions in W31, the kinematic analysis in section 3.3 



also indicative of evolutionary differences between the two com- 
plexes in general. The virial mass as presented in section [33] is 
directly proportional to the observed line width (eq. [T]). In the 
simple virial equilibrium picture, having a ratio Myk/Mpeak ~ 1 
would imply that the clumps could be stable against collapse. 
However, except for a few notable sources like G10.3W2 (Fig.|3] 
and clump 6 in Fig. |4] see also Table |2] and Fig.[TT]) which are 
candidates for starless clumps, most other regions in this com- 
plex show various signs of star formation. Hence, they are un- 
likely to be in virial equilibrium. It is rather the opposite, the 
feedback from ongoing star formation (e.g., outflows or radia- 
tion) significantly broadens the observed line width and the virial 
assumption is not properly applicable under such circumstances. 



While we do not know whether the clumps with high 
Mvir/Mpeak ratios are still bound or already expanding again from 
the inner energy sources, the fact that we see ongoing star forma- 
tion indicates that they are still at least partly bound. Large virial 
parameters are also consistent with pressure-confined clumps 
as discussed by|Bertoldi & McKee| (|1992| ) and more recently 
by |Lada et al.| ( |2QQ8| ), [Dobbs et al.| pOTTf and Kainulainen et 
al. (subm.). High virial mass ratios have also been found by 
[Simon et al.| r20Ql ) for fo ur molecular cloud s studied via the 
Galactic Ring Survey (GRS, Jackson et al. 2006 ). It is interesting 
to note that the three less evolved clouds of their sample have on 
average even higher virial mass ratios than their most luminous 
source, the galactic mini-starburst W49. The peak of their ratio 
between virial mass and LTE mass (derived via LTE calculations 
for the ^^CO(l-O) emission) for W49 is around 2, whereas we 
find an average of Myir/Mpeak for W31 of ~ 1.1. Hence, similar 
to W49 studied by |Simon et al.[ ( |2001| , the clumps in W31 are 
in the regime of being bound. 

Maybe more surprising is that all high-mass clumps associ- 
ated with GlO.3-0.1 have Myir/Mpeak ratios below 1. While this 
could be expected for the starless clump candidate G10.3W2 if 
it were in an unstable state just at the verge of collapse and had 
no additional turbulent support, it is less expected for the other 
active star- forming clumps where multiple additional line width 
broadening mechanisms should be at work, e.g., outflows and/or 
radiation. Nevertheless, low Myir/Mpeak ratios clearly imply a 
comparatively small turbulent line width with respect to most 
of the clumps associated with the southern region Gl 0.2-0.3. 
Furthermore, magnetic fields may be more important in these 
supposedly younger clumps acting as additional source of sta- 
bility, also implying low Myir/Mpeak ratios. 

While the overall characteristics of the two complexes 
GlO.2-0.3 and GlO.3-0.1 appear similar - both have similar lu- 
minosities (see section[T]), strong submm and mid-infrared emis- 
sion as well as an associated near-infrared cluster - there is 
an apparent evolutionary difference between the two. In addi- 
tion to the above discussed different Myir/Mpeak ratios between 
GlO.2-0.3 and GlO.3-0.1, the class O/I and class II sources 
around the southern GlO.2-0.3 region appear considerably more 
widely distributed than those around the northern region G10.3- 
0.1 (see section [3JrT] ). Furthermore, it is also interesting to note 
that only GIO. 3-0.1 exhibits several sites of Cla ss ii CH3QH 
masers whereas none is found toward GIO. 2-0.3 ( [Walsh et al.| 
1998). Since these CH3OH masers are also usually found to- 
ward younger regions (e.g.. Fish 2007), the presence of them 
toward only the northern region GIO. 3-0.1 combined with the 
above discussed additional differences strongly suggests that the 
northern GIO. 3-0.1 complex on average should be younger than 
the southern G 10.2-0. 3 region. The latter then had more time al- 
ready to stir up the gas of the environment and hence increase 
the observed line width. 

On top of all these evolutionary effects, we also find a closer 
association of the IRAC-identified younger class O/I sources 
with the central high-mass gas clumps while the more evolved 
class II sources appear more broadly distributed over the whole 
field (Fig.[T]). 

4.2. The shape of the clump mass function 

The power-law index a ~ 1.5 ± 0.3 of the clump mass func- 
tion discussed in Fig. [5] (AA/^/ AM oc M"^) is consistent with the 
clump mass function derived for molecular clouds by molec- 
ular line CO observations (a ~ 1.6, e.g., [Stutzki & Guesten 
[T99Q1 [BTite|[T993l [Kramer et al.l[T998l [Simon et al.[[2001D , but 
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it is considerably flatter than distributions derived for other star- 
forming regions from dust continuum observations or extinction 
mapping, which more resemble the Salpeter slope (a - 2.35, 



e.g., 


Motte et al.| 


1998; Beltran et al. 2006; Reid & Wilson 


2006 


Alves et al.| 


2007). In the past, the diff'erences in power- 



law mass distributions derived from the diff'erent tracers (CO 
versus dust emission/extinction) was often attributed to the dif- 
ferent density regimes they are tracing: The CO observations 
are more sensitive to the diff'use and transient gas whereas 
the dust emission/extinction rather traces the denser, gravita- 
tionally or pressure-bound cores. In this picture, the diff'erent 
power-law distributions would reflect a structural change of the 
cloud/clump properties from transient to bound structures. 

However, this picture does not hold for our data of the W3 1 
complex since we know that the region is actively star-forming 
and not a transient structure. Since cluster mass functions are 
also flatter than stellar mass functions (power-law distributions 
of 2.0, e.g., Bik et al.|2003l or |de Grijs et al.|2003| versus 2.35 of 
the Salpeter slope, respectively), one possibility to explain the 
discrepancies is that we are tracing high-mass cluster-forming 
regions at far dista nce s like W31 wherea s the observations by 
[Motte et aL] ( |1998| ) or [Alves et al.| ( |2007| ) dissect low-mass re- 
gions like p Ophiuchus or the Pipe at distances closer < 150 pc. 
Unfortunately, this explanation is not consistent with all obser- 
vations since the mass distributions derived by Beltran et al. 
( [2006[ ) or [Reid & Wi lson ( 2006) are also targeting high-mass 
star-forming regions at distances of several kpc. 

Another possibility for a flatter clump mass function with re- 
spect to the IMF could be that the star formation efficiency (SFE) 
varies with clump mass. Higher mass gas clumps may exhibit on 
average a lower SFE which would steepen the slope from the 
clump mass function to the I MF. Similar results were recently 
inferred byjParmentier (20111 In a similar direction, it is worth 
noting that [Simon et al.[ |2()01) also inferred for their most lumi- 
nous mini-starburst W49 the flattest clump mass function with a 
power-law index of 1.56, in very close agreement with our result 
for the high-luminosity region W3 1 . Therefore, there are several 
indications that the most luminous regions could also exhibit the 
flattest clump mass functions. 

One potentially more technical explanation for the relatively 
flat clump mass function we derive for W3 1 could be that our 
used uniform temperature of 50 K more severely afl'ects a than 
we anticipated. Higher temperatures decrease the mass estimates 
of a clump whereas lower temperatures raise the mass esti- 
mates. In this picture, it could well be that the higher mass 
clumps closer to the centers of the region have higher temper- 
atures whereas lower mass clumps may still be colder. This ef- 
fect would steepen the mass function. Although with the current 
data we cannot properly quantify this efl'ect, it is unlikely that 
the mass function would steepen to a Salpeter slope. 

Therefore, we suggest that the observed flat mass function 
reflects the facts that the resulting cluster mass functions are also 
flatter than the stellar initial mass function, and furthermore that 
the SFE may vary with clump mass. 

5. Conclusions 

The most luminous and high-mass star formation appears to take 
place in a highly structured and non-uniform fashion. The multi- 
wavelength continuum and spectral line analysis of the difl'er- 
ent sub-regions within the very luminous W31 regions reveals 
evolutionary difl'erences on large scales between the two sub- 
complexes GlO.2-0.3 and GlO.3-0.1, but also within each of 
these regions on much smaller scales. While many clumps have 



virial mass ratios close to unity, we also find low- turbulence, po- 
tentially still starless high-mass gas clumps that can reside in 
almost the direct vicinity of ongoing and already finished high- 
mass cluster-forming regions. Furthermore, these dense active 
gas clumps appear to be surrounded by a halo-like distribution 
of already more evolved class 0/1 and class II sources. We also 
find a tight correlation between the warm dust tracing 24/dm 
emission and the ionized gas tracing cm emission, implying that 
warm dust at temperatures around 1000 K can spatially coexist 
with ionized gas in the 10"^ K regime. In addition to this, our data 
indicate that the clump mass function of W3 1 is considerably 
flatter (power-law index a ~ 1.5 + 0.3) than the initial mass func- 
tion, but it is consistent with the mass functions derived for other 
molecular clouds. Since these gas clumps trace cluster scales, 
this is consistent with the flatter cluster mass function compared 
to the stellar mass function. Furthermore, the analysis tentatively 
suggests that the star formation efficiency may decrease with in- 
creasing clump mass. 
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Table 1. Clump parameters for original data. The first part shows the coordinates, and the second and third parts the SlS/im and the 
C^^0(2-l) data, respectively. 



# 


R.A. 


Dec. 


»^peak 


S int 


^efF 


(dust) 


M 


^peak 




Vpeak 


Av 


A^H2 

(C^^O) 






J2000.0 


J2000.0 


Jy 


Jy 


pc 




Mo 


K 


Kkms-i 


kms~^ 


kms~^ 




Mo 


1 


18 08 55.83 


-20 05 56.1 


7.16 


27.37 


1.2 


34.3 


7359 


2.8 


7.2 


11.2 


2.4 


2.8 


285 


2 


18 08 59.96 


-20 03 36.4 


4.93 


16.86 


1.2 


23.6 


4534 


1.3 


2.6 


10.1 


1.8 


1.0 


167 


3 


18 09 01.61 


-20 05 09.6 


4.36 


30.45 


1.7 


20.9 


8187 


2.1 


6.5 


11.9 


2.9 


2.5 


416 


4 


18 09 24.34 


-20 15 38.5 


4.31 


21.31 


1.3 


20.6 


5731 














5 


18 08 49.22 


-20 05 56.1 


4.02 


11.41 


0.8 


19.2 


3067 


0.8 


4.2 


11.7 


4.8 


1.6 


1175 


6 


18 08 46.33 


-20 05 50.2 


3.81 


18.61 


1.3 


18.2 


5003 


1.7 


1.9 


10.6 


1.0 


0.7 


50 


7 


18 09 27.65 


-20 19 08.1 


3.68 


27.65 


1.3 


17.6 


7434 


1.8 


9.3 


10.1 


4.8 


3.6 


1159 


8 


18 09 21.45 


-20 19 31.4 


3.43 


25.94 


1.5 


16.4 


6974 


2.5 


11.7 


10.2 


4.3 


4.5 


938 


9 


18 08 52.11 


-20 06 07.8 


2.37 


10.44 


1.1 


11.3 


2808 


2.2 


6.5 


11.1 


2.8 


2.5 


405 


10 


18 09 26.82 


-20 17 29.1 


2.33 


15.49 


1.3 


11.1 


4164 


1.0 


6.5 


11.9 


6.3 


2.5 


2014 


11 


18 09 25.58 


-20 18 27.3 


2.31 


17.87 


1.3 


11.1 


4804 


1.0 


7.6 


12.9 


6.9 


2.9 


2377 


12 


18 09 28.89 


-20 16 48.3 


2.30 


18.61 


1.4 


11.0 


5004 


1.0 


6.7 


11.9 


6.2 


2.6 


1909 


13 


18 09 23.09 


-20 08 04.3 


2.28 


8.06 


1.2 


10.9 


2168 














14 


18 09 03.26 


-20 03 01.5 


2.09 


13.17 


1.3 


10.0 


3540 


0.8 


1.2 


11.2 


1.4 


0.5 


101 


15 


18 09 00.36 


-20 11 33.9 


1.88 


16.27 


1.7 


9.0 


4374 














16 


18 09 20.62 


-20 15 03.6 


1.65 


3.93 


0.9 


7.9 


1055 














17 


18 09 35.92 


-20 18 44.7 


1.53 


6.19 


1.0 


7.3 


1665 


0.8 


2.6 


9.7 


3.1 


1.0 


469 


18 


18 09 33.85 


-20 17 52.4 


1.46 


6.92 


1.1 


7.0 


1860 














19 


18 08 47.56 


-20 06 48.5 


1.44 


7.31 


1.1 


6.9 


1967 


1.0 


3.4 


9.7 


3.1 


1.3 


499 


20 


18 09 35.56 


-20 21 27.8 


1.30 


3.09 


0.7 


6.2 


830 














21 


18 09 21.03 


-20 16 07.6 


1.26 


11.90 


1.4 


6.0 


3201 














22 


18 09 26.83 


-20 21 22.0 


1.22 


4.29 


0.9 


5.8 


1154 


1.8 


7.1 


11.1 


3.7 


2.7 


687 


23 


18 08 57.89 


-20 07 06.0 


1.08 


5.35 


1.1 


5.2 


1439 


1.5 


6.6 


11.3 


4.3 


2.5 


923 


24 


18 09 24.76 


-20 20 47.1 


1.03 


5.53 


1.0 


4.9 


1486 


3.0 


15.6 


11.4 


4.8 


6.0 


1173 


25 


18 09 23.09 


-20 09 31.6 


1.03 


4.39 


0.9 


4.9 


1180 














26 


18 09 14.42 


-20 18 56.5 


0.99 


5.32 


1.1 


4.7 


1431 


5.4 


21.6 


12.7 


3.7 


8.3 


702 


27 


18 08 41.36 


-20 07 35.0 


0.90 


3.26 


0.9 


4.3 


862 


1.5 


4.2 


9.1 


2.7 


1.6 


368 


28 


18 09 06.57 


-20 03 18.9 


0.89 


4.22 


1.1 


4.3 


1134 


1.1 


2.0 


11.2 


1.7 


0.8 


153 


29 


18 09 05.33 


-20 03 42.2 


0.89 


9.02 


1.4 


4.3 


2425 


1.4 


2.5 


11.0 


1.7 


0.9 


137 


30 


18 09 29.72 


-20 20 06.3 


0.87 


2.82 


0.8 


4.2 


758 


0.9 


8.2 


10.3 


8.1 


3.1 


3331 


31 


18 08 39.26 


-20 18 56.3 


0.83 


2.06 


0.8 


4.0 


553 














32 


18 09 21.02 


-20 02 03.2 


0.81 


7.00 


1.4 


3.9 


1881 














33 


18 09 39.23 


-20 19 31.3 


0.80 


2.05 


0.7 


3.8 


552 


0.7 


0.8 


6.6 


1.0 


0.3 


50 


34 


18 09 12.76 


-20 17 46.6 


0.77 


2.09 


0.8 


3.7 


563 


1.3 


11.0 


12.5 


7.8 


4.2 


3099 


35 


18 09 34.69 


-20 22 08.6 


0.72 


1.25 


0.6 


3.4 


337 














36 


18 09 37.16 


-20 16 42.4 


0.71 


4.90 


1.2 


3.4 


1318 














37 


18 09 25.18 


-20 25 44.1 


0.70. 


4.88 


1.3 


3.3 


1312 


2.5 


11.9 


11.2 


4.6 


4.6 


1055 


38 


18 09 17.30 


-20 07 11.9 


0.69 


2.80 


0.9 


3.3 


752 


0.9 


2.6 


13.4 


2.6 


1.0 


347 


39 


18 09 35.51 


-20 20 52.9 


0.68 


1.75 


0.7 


3.3 


471 


0.8 


3.0 


10.5 


3.3 


1.1 


559 


40 


18 09 34.28 


-20 22 26.0 


0.67 


1.02 


0.5 


3.2 


273 


0.5 


1.4 


9.8 


2.6 


0.6 


331 


41 


18 08 59.95 


-20 06 31.1 


0.66 


2.52 


0.8 


3.2 


677 


1.7 


2.3 


11.7 


1.3 


0.8 


79 


42 


18 09 13.17 


-20 18 04.1 


0.66 


1.85 


0.7 


3.2 


498 


1.7 


10.4 


13.2 


5.8 


4.0 


1722 


43 


18 09 31.80 


-20 23 59.2 


0.66 


2.35 


0.8 


3.2 


632 


0.9 


2.4 


16.1 


2.5 


0.9 


313 


44 


18 09 26.40 


-20 14 28.6 


0.65 


3.57 


1.0 


3.1 


961 


0.5 


1.3 


9.8 


2.2 


0.5 


255 


45 


18 09 37.99 


-20 19 48.8 


0.65 


0.76 


0.5 


3.1 


214 














46 


18 09 21.85 


-20 10 00.7 


0.64 


1.72 


0.7 


3.1 


461 














47 


18 09 20.21 


-20 20 29.6 


0.64 


3.29 


1.0 


3.1 


885 


6.3 


25.9 


12.2 


3.9 


10. 


747 


48 


18 08 55.40 


-20 1104.8 


0.60 


4.84 


1.3 


2.9 


1301 


1.2 


2.3 


14.0 


1.8 


0.9 


158 


49 


18 08 36.36 


-20 19 37.0 


0.57 


2.23 


0.9 


2.7 


599 














50 


18 09 14.83 


-20 15 32.7 


0.56 


1.97 


0.8 


2.7 


529 














51 


18 09 33.02 


-20 13 59.4 


0.56 


1.29 


0.6 


2.7 


348 


0.9 


2.3 


9.0 


2.6 


0.9 


331 


52 


18 09 01.58 


-20 30 23.6 


0.56 


3.43 


1.1 


2.7 


923 


0.9 


1.4 


11.7 


1.4 


0.5 


100 


53 


18 08 45.08 


-20 07 35.0 


0.54 


3.05 


1.0 


2.6 


820 














54 


18 09 09.86 


-20 29 37.1 


0.53 


1.43 


0.7 


2.5 


385 


0.8 


2.9 


10.7 


3.4 


1.1 


568 


55 


18 08 41.36 


-20 08 21.6 


0.51 


1.58 


0.7 


2.4 


426 


0.8 


2.2 


10.4 


2.7 


0.9 


355 


56 


18 09 06.97 


-20 26 36.5 


0.50 


1.76 


0.8 


2.4 


474 


0.8 


0.9 


11.7 


1.1 


0.4 


66 


57 


18 08 31.40 


-20 18 56.2 


0.50 


1.54 


0.7 


2.4 


414 














58 


18 09 15.24 


-20 16 36.7 


0.48 


1.65 


0.7 


2.3 


444 


1.0 


4.0 


9.7 


3.9 


1.5 


762 


59 


18 09 11.52 


-20 16 19.3 


0.47 


1.26 


0.7 


2.2 


339 














60 


18 09 37.17 


-20 20 23.7 


0.47 


0.57 


0.4 


2.2 


154 














61 


18 08 40.54 


-20 08 21.6 


0.46 


0.62 


0.5 


2.2 


168 


1.2 


3.6 


10.0 


2.9 


1.4 


410 


62 


18 08 38.01 


-20 20 58.6 


0.45 


0.88 


0.6 


2.2 


231 














63 


18 09 09.04 


-20 25 55.8 


0.45 


1.63 


0.8 


2.2 


439 














64 


18 09 13.59 


-20 16 19.3 


0.45 


0.88 


0.6 


2.2 


238 














65 


18 09 41.28 


-20 12 08.7 


0.44 


1.97 


0.9 


2.1 


529 














66 


18 08 40.53 


-20 08 33.2 


0.43 


1.08 


0.6 


2.1 


289 














67 


18 08 45.91 


-20 07 05.9 


0.43 


1.17 


0.6 


2.1 


314 


0.8 


1.4 


8.8 


1.7 


0.5 


152 


68 


18 08 43.02 


-20 06 54.3 


0.43 


2.49 


0.9 


2.1 


669 


1.1 


3.6 


15.8 


3.2 


1.4 


505 


69 


18 09 15.66 


-20 16 19.3 


0.42 


0.76 


0.5 


2.0 


203 














70 


18 08 43.87 


-19 58 56.8 


0.41 


0.93 


0.6 


2.0 


249 














71 


18 09 06.55 


-20 28 38.8 


0.39 


0.77 


0.5 


1.9 


206 


1.0 


2.9 


10.6 


2.9 


1.1 


413 


72 


18 08 21.87 


-20 22 54.8 


0.32 


0.45 


0.4 


1.5 


121 
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Table 2. Clump parameters for 27.5'' resolution. The first part shows the coordinates, and the second and third parts the SlSjim and 
the C^^0(2-l) data, respectively. 875 fim clump parameters for 27.5'' resolution. 



# 


R.A. 


Dec. 


»^peak 


S int 


^efF 


(dust) 


M 


Mpeak 


^peak 


J Tp^akdv 


Vpeak 


Av 


A^H2 






J2000.0 


J2000.0 


Jy 


Jy 


pc 




Mo 


Mo 


K 


Kkms-1 


kms~^ 


kms~^ 




Mo 


1 


18 08 55.83 


-20 06 02.0 


11.22 


36.16 


1.7 


26.2 


9723 


3017 


3.0 


7.5 


11.3 


2.4 


2.9 


279 


2 


18 08 59.96 


-20 03 36.4 


6.75 


16.31 


1.2 


15.7 


4386 


1815 


1.5 


2.6 


10.0 


1.7 


0.9 


141 


3 


18 09 24.34 


-20 15 38.5 


6.71 


31.44 


2.0 


15.6 


8454 


1804 














4 


18 09 01.61 


-20 05 09.6 


6.67 


34.26 


2.2 


15.6 


9213 


1793 


1.8 


6.4 


11.7 


3.3 


2.5 


538 


5 


18 09 27.65 


-20 19 08.1 


6.50 


36.43 


1.7 


15.2 


9796 


1749 


1.8 


9.0 


10.0 


4.7 


3.5 


1129 


6 


18 08 46.74 


-20 05 50.2 


6.44 


38.99 


2.1 


15.0 


1048 


1732 


1.1 


2.6 


10.5 


2.2 


1.0 


244 


7 


18 09 21.45 


-20 19 31.4 


5.68 


33.06 


2.1 


13.3 


8890 


1528 


2.4 


11.8 


10.4 


4.6 


4.5 


1052 


8 


18 09 26.41 


-20 17 40.8 


4.27 


23.00 


1.7 


10.0 


6183 


1147 


1.3 


7.4 


12.4 


5.6 


2.9 


1568 


9 


18 09 28.89 


-20 16 48.3 


3.77 


21.92 


1.8 


8.8 


5894 


1015 


1.5 


8.1 


12.0 


5.1 


3.1 


1336 


10 


18 09 03.26 


-20 03 07.3 


3.59 


23.71 


2.1 


8.4 


6376 


965 


0.6 


0.9 


11.3 


1.4 


0.3 


99 


11 


18 09 00.36 


-20 11 33.9 


3.22 


16.32 


1.9 


7.5 


4388 


866 














12 


18 09 23.09 


-20 08 04.3 


3.10 


7.36 


1.3 


7.2 


1979 


833 














13 


18 09 36.34 


-20 18 44.7 


2.44 


6.03 


1.1 


5.7 


1620 


656 


0.6 


2.9 


9.3 


4.6 


1.1 


1074 


14 


18 09 33.85 


-20 17 52.4 


2.34 


11.23 


1.8 


5.5 


3021 


629 














15 


18 09 20.62 


-20 15 03.6 


2.17 


11.43 


1.8 


5.1 


3074 


585 


0.6 


1.5 


10.3 


2.3 


0.6 


273 


16 


18 09 26.83 


-20 21 22.0 


1.95 


4.11 


1.0 


4.5 


1105 


524 


2.1 


1.1 


11.2 


3.5 


3.0 


602 


17 


18 09 35.52 


-20 21 27.8 


1.89 


7.37 


1.5 


4.4 


1980 


507 














18 


18 08 58.30 


-20 07 06.0 


1.81 


7.10 


1.4 


4.2 


1908 


485 


1.2 


5.5 


11.3 


4.4 


2.1 


959 


19 


18 09 25.17 


-20 20 47.1 


1.74 


7.48 


1.4 


4.1 


2011 


469 


3.2 


16.2 


11.5 


4.8 


6.2 


1169 


20 


18 09 24.33 


-20 09 08.3 


1.62 


6.94 


1.4 


3.8 


1865 


436 














21 


18 08 41.36 


-20 07 35.0 


1.42 


12.96 


2.1 


3.3 


3483 


381 


1.2 


3.6 


9.2 


2.8 


1.4 


386 


22 


18 09 21.02 


-20 02 03.2 


1.33 


7.63 


1.7 


3.1 


2052 


359 














23 


18 09 39.23 


-20 19 31.3 


1.15 


3.38 


1.1 


2.7 


910 


309 


0.6 


0.9 


6.7 


1.5 


3.5 


108 


24 


18 09 13.17 


-20 17 52.4 


1.09 


10.21 


2.0 


2.5 


2745 


292 


1.2 


2.9 


10.4 


2.3 


1.1 


265 


25 


18 09 17.30 


-20 07 11.9 


1.09 


2.70 


1.0 


2.5 


727 


292 


1.3 


2.7 


13.5 


2.0 


1.0 


208 


26 


18 08 38.85 


-20 18 56.3 


1.07 


2.09 


0.9 


2.5 


562 


287 














27 


18 09 32.21 


-20 23 59.2 


1.03 


2.88 


1.1 


2.4 


775 


276 


0.9 


2.7 


15.9 


2.8 


1.1 


385 


28 


18 09 25.18 


-20 25 44.1 


1.01 


4.63 


1.5 


2.3 


1244 


270 


2.6 


12.1 


11.1 


4.3 


4.7 


939 


29 


18 08 54.99 


-20 1104.8 


0.98 


5.04 


1.4 


2.3 


1356 


265 


2.6 


14.3 


2.3 


1.0 


262 




30 


18 09 01.58 


-20 30 29.4 


0.94 


3.57 


1.2 


2.2 


960 


253 














31 


18 09 32.60 


-20 14 05.3 


0.86 


1.80 


0.9 


2.0 


484 


232 


0.7 


1.7 


8.3 


2.2 


0.7 


247 


32 


18 09 07.38 


-20 26 30.7 


0.84 


5.53 


1.6 


2.0 


1486 


226 














33 


18 08 36.36 


-20 19 31.2 


0.80 


2.38 


1.0 


1.9 


641 


215 














34 


18 08 32.23 


-20 19 13.7 


0.78 


2.01 


1.0 


1.8 


540 


210 














35 


18 09 09.45 


-20 29 37.1 


0.78 


1.69 


0.9 


1.8 


456 


210 


1.0 


2.9 


10.8 


2.8 


1.1 


398 


36 


18 09 28.49 


-20 25 15.0 


0.68 


3.44 


1.3 


1.6 


924 


182 














37 


18 09 40.87 


-20 12 26.2 


0.66 


2.42 


1.1 


1.5 


652 


177 














38 


18 08 38.01 


-20 20 58.6 


0.64 


1.34 


0.9 


1.5 


361 


171 














39 


18 08 43.46 


-19 58 50.9 


0.59 


1.31 


0.9 


1.4 


352 


160 














40 


18 08 21.87 


-20 22 54.8 


0.51 


0.60 


0.6 


1.2 


160 


138 














41 


18 08 15.71 


-20 13 47.3 


0.43 


0.50 


0.6 


1.0 


133 


116 














42 


18 08 03.82 


-20 01 33.3 


0.41 


0.34 


0.5 


1.0 


91 


110 














43 


18 08 31.04 


-20 05 44.2 


0.39 


0.50 


0.6 


0.9 


135 


105 














44 


18 09 38.35 


-20 00 53.2 


0.39 


0.35 


0.5 


0.9 


93 


105 















The table parameters are 875yLim peak and integrated fluxes 5 peak and 5 int, eff'ective linear clump radius from the clumpfind search reff, and 
calculated H2 column densities A^H2 and gas masses M and Mpeak, derived from the total flux and peak flux, respectively. The C^^0(2-l) part 
shows the peak and integrated line intensities Tpeak and J Tpeak^^v, the peak velocities Vpeak and line widths Av, and the derived H2 column densities 
A^H2 and virial gas masses Myir. 



